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SUMMARY 

Combustion time lag is treated in a form that 
includes its dependency on injection velocity , such 
dependency being indicated by recent studies . This 
generalized form is used in analyses of low-frequency 
chamber dynamics and combustion instability (tank- 
fed rocket). The theoretical responses and stability 
boundaries obtained are compared with those given 
with the previous time-lag concept . It is concluded 
that the injection-velocity effect on the time lag 
cannot be neglected in the theory of chamber dynamics 
and combustion instability 

INTRODUCTION 

In the field of liquid -propellant rocket engines, 
there is widespread interest in the theory of 
combustion-chamber dynamics. In the past, this 
interest was generated primarily by the frequent 
occurrence of combustion instability. Often 
costly and time-consuming development programs 
wore required for new designs in order to over- 
come this destructive phenomenon. At the pres- 
ent time, although combustion instability is still 
of much concern, increasing attention is being 
given to defining the chamber dynamics for control 
purposes. This late emphasis on controls results 
from the stringent trajectory requirements of 
space missions. 

The foundation of existing theory was laid in 
1941 when Von Ktfrmdn recognized that a time 
lag between propellant injection and combustion 
was instrumental in combustion instability (ref. 
1). This combustion time lag (also known as 
dead time) was first treated as being of constant 
duration. It was used in this manner in analyses 
by Gunder and Friant (ref. 2), Yachter and 
Waldingcr (ref. 3), Summerfield (ref. 4), and Tee, 


Gore, and Ross (ref. 5). Crocco (refs. G and 7) 
improved the concept by considering the time lag 
to he a variable with a duration dependent on 
changing conditions along the path of the pro- 
pellant in the chamber. For simplicity he corre- 
lated the various conditions to chamber pressure. 
Crocco’s time-lag concept was used in stability 
studies by Tsien (ref. 8) and Marble and Cox 
(ref. 9). At the present time, Crocco, Grey, and 
Matthews (ref. 10) are using this concept of 
pressure-dependency to determine time lags from 
experimental frequency response. This same con- 
cept, is the basis of the extensive theory of com- 
bustion instability developed by Crocco and 
Cheng (ref. 11). 

Recent developments, however, indicate the 
time lag must be treated as a function of injection 
velocity as well as conditions after injection. 
Studies have shown that the atomization of t lie 
propellant is sensitive to injection velocity (refs. 
12 and 13). Propellant atomization, in turn, lias 
long been considered an important factor in 
determining the duration of the time lag (ref. 11, 
e.g.). Substantiation of this belief is given by 
the correlation of engine performance with vapori- 
zation in reference 14, which indicates that the 
atomization-sensitive vaporization process may 
determine the time lag. Experimental evidence 
that directly shows the variation of time lag with 
injection velocity is reported by Penner and Fuhs 
(ref. 15). 

The time-lag concept, therefore, is generalized 
in the present treatment to include the effect of 
injection velocity. This is done by treating the 
time lag as being dependent on initial propellant 
atomization as well as conditions encountered by 
the propellant in the chamber; the initial atomi- 
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zation is considered to be a function of injection 
velocity. The chamber variables along the path 
of the propellant after injection are correlated, in 
the manner of CTocco, to chamber pressure. 

Analyses are presented that show the importance 
of using the more general time-lag concept in the 
treatment of chamber dynamics and combustion 
instability. With the generalized concept, a 
small-perturbation equation of the combustion 
chamber is derived for a bipropellant rocket. The 
equation treats the range of low frequencies in- 
volved in the combustion instability called chug- 
ging. This is also the frequency range of primary 
interest in controls work. From the chamber 
equation, the transfer functions are formulated, 
and the frequency response is presented for pos- 
sible sensitivities of the time lag to injection 
velocity. The chamber equation is then used in 
analyzing the stability of a tank-fed rocket in 
order to illustrate the injection-velocity role in 
combustion instability. 

SYMBOLS 

/? mass-transfer number 

D drop diameter 

E initial or final level of time-lag processes 

FJ rate of time-lag processes 

h inject ion -hole diameter, in. 

K pl2(P oJ -p) 

mass of burnt gas in combustion chamber 
m exponent in F ~ D m p n 

n exponent in F^D m p n (interaction index) 

P of tank pressure of oxidizer or fuel 

p pressure of gas in combustion chamber 

Q injection-velocity sensitivity, 

{ zero (eq. (20)) 

(r+l)(3— m)ya (eq. (21)) 

(3 —m)ya (eqs. ((22) and (25)) 
R gas constant 

Re Reynolds number 

r mixture ratio, wjw f 

Sc Schmidt number 

s dimensionless complex operator 

T gas temperature 

t time 

U drop velocity relative to gas in chamber 

V injection velocity 

ic b gas generation rate 

w e exhaust rate 

w f _ fuel rate 

w t total injection rate 


w 0 oxidizer rate 

■ '/{'-Gfen 

-ii'-mry 

z dimensionless lime, t/6 g 

a atomization index 

ft dimensionless frequency of neutral oscilla- 

tion, d K w } where o> is angular frequency 
A( ) small perturbations 
d g gas residence time 

fi gas viscosity 

pi liquid density 

a time lag (dead time) 

r dimensionless time lag, ajd g 

Superscript: 

( ) steady -state values 

COMBUSTION-CHAMBER DYNAMICS 
THE TIME LAG 

The combustion in a liquid-propellant rocket 
engine is very complicated and still not fully 
understood. Intricate processes such as atomi- 
zation, mixing, vaporization, and chemical re- 
action take place before the final products of 
combustion evolve from the injected liquids. 
Although the evolution is certainly gradual in 
both time and space, a treatment of the phe- 
nomenon in its gradual nature is not warranted 
by the present knowledge of the processes involved 
and of the chamber conditions affecting these 
processes. A working model is obtained, there- 
fore, in the manner of previous authors by con- 
sidering small fractions of the final products to be 
generated in a discontinuous manner. That is, 
the gradual evolution of small elements of pro- 
pellants into a particle of burnt gas is replaced 
by a sudden conversion. The time interval be- 
tween injection and the sudden conversion is 
called time lag. In general, different particles of 
burnt gas will have different time lags, and the 
time lag of each will be a time-dependent quantity. 
Nevertheless, considerable simplification is 
achieved in that the complicated combustion 
process is represented merely by a distribution of 
gas sources. 

Tt has been shown (ref. 11) that the rate at 
which the burnt gas is generated from these 
sources is dependent on the rate of change of the 
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associated time lags. It is necessary, therefore, to 
consider the manner in which the time lag varies 
before proceeding with the analysis of the chamber 
dynamics. 

Previous authors have assumed that the vari- 
ations in the time lag are caused only by variations 
in the physical factors encountered by the pro- 
pellant after injection into the chamber. The as- 
sumption is tacitly contained in the equation 

f f(t')dt' = E(t) = const. (1) 

J 1-0 

used to define the dependency of the time lag cr on 
the physical factors. The integral, which is taken 
over the path of the propellant element from the 
time of injection t'—t—a to the time of burning 
t'—t, represents the accumulation during the time 
lag of the processes determining its duration. The 
process rate f(t r ) is influenced by the chamber 
variables such as pressure, temperature, and rela- 
tive velocities between propellant and gas. The 
constant right-hand member E(t) represents the 
final level of the processes, that is, the level to 
which they must accumulate before the sudden 
conversion into burnt gas. The equation, of 
course, is sufficient only with the assumption that 
the initial level is constant, the general form being 


one having longer drop lifetimes.) The rate at 
which the drop size is reduced at each instant 
along its path can be expressed by the pro- 
portionality 

dZ> 3 

(m<3) 

where F represents a function dependent on 
chamber variables. The relation can also be 
written as 


dD 3-w 

dt' 


(3) 


with the proportionality constant being absorbed 
in the new function/. By specifying the size of 
the drop at the end of the time lag as D(t ), the 
duration of the time lag is defined from equation 
(3) as follows: 

which, upon integration of the left-hand member, 
becomes 



f(t')dt'=IP- m {t) 


(5) 


E(t— <t)+ f /(f')d£' = 7?(0 = const. (2) 
J t-a 

where the term E(t a) represents the level of the 
processes, or the condition of the propellant, at 
the instant of injection. 

In the following analysis, this assumption is 
dropped, and the condition of the propellant at 
the instant of injection is treated as a variable 
with an effect on the duration of the time lag. To 
do this quantitatively, of course, requires that the 
initial (and final) condition in equation (2) be 
defined explicitly. 

The more explicit equation is obtained in the 
following manner. The liquid drop is treated as 
the elemental unit of propellant, and it is postu- 
lated that its actual combustion can be adequately 
depicted by a time-lag approximation controlled 
by the lifetime of the drop, or a portion thereof. 
In other words, the sudden conversion to burnt 
gas, the end of the time lag, is considered to occur 
when the drop has been reduced to a definite size, 
(For bipropellants, the species to be treated is the 


It should be noted that the treatment result- 
ing in equation (5) does not necessarily imply 
the actual combustion process to be one of “drop- 
burning, 55 in which the vapors issuing from the 
drop are immediately burned. The model suf- 
fices, rather, if the vapors formed along the path 
of the drop accumulate and bum at a later time; 
the only condition being, as previously stated, 
that the actual combustion can be adequately 
approximated by a discontinuous transformation 
taken at a definite drop size. The model may 
not be sufficiently general, of course, for some 
applications. More generality would be obtained 
by treating the time lag as the sum of the drop 
lifetime and a gas-phase delay, each with its own 
process rate. For high-frequency applications, in 
which even small delays are significant, the more 
general treatment appears essential. In fact, the 
recent work of Crocco, Grey, and Ilarrje (ref. 16) 
shows conclusively that small delays highly de- 
pendent on chamber pressure are instrumental 
in high-frequency combustion instability (scream- 
ing). Such delays are probably in the gas phase. 
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For low-frequency applications, such as those 
in the later sections of the present report, the 
general formulation may be required if there is 
considerable 1 - combustion delay after the drop 
vanishes. It is not deemed appropriate, however, 
to complicate the present work by this generaliza- 
tion since interest is presumed to be centered 
primarily in rockets with good liquid phase mixing. 

In proceeding, therefore, the rate of change of 
time lag is found from equation (5) to be 


<!*_ i (f > m-iit-c) 

df jit — <t) df 1 j j{t-c) 


(G) 


which, for the small-port nrhution treatment, can 
be written as 


da 
d t 


i dzy— 

fit -a) d t 


«- 5 ) 


4/(0 — A f(t-a) Kt)-j{t-a) 
Jit — a) J{t—a) 


(7) 


The process rale /can be expressed, in the manner 
of Croc co, as 

f^P" 


where n } the interaction index, includes the 1 effects 
of other chamber variables under the assumption 
they can be correlated to pressure. The rate 
perturbation, therefore, is related to the pressure 
perturbation by 

A / Ap 

7~ n P 

and, with uniform pressure in steady-stale, 
equation (7) becomes 


dcr 1 (]D 3 ~ m Ap(f)-Ap(f-Z) 

di = -j—\r — j 


(8) 


An expression for the rale / can be obtained by 
integrating equation (5) for steady-state condi- 
tions; substitution into equation (S) yields 

o_ (ID 3 w 

dt~~ D*- m {t-a)-D*- m {t) d t ' ~ a) 

Agffl-AML =D (9) 

P 

The dependence of initial drop size on injec- 
tion velocity indicated by references 12 and 13 


can be expressed for the small perturbation 
treatment as 

D{1— <r) ^V~ a (t — <r) 

Using this expression with equation (9) gives the 
following equation relating the rate of change of 
time lag to perturbations in injection velocity 
and chamber pressure: 



It should be rioted that the pressure perturbation 
Ap(i) appearing in this equation is a function of 
time and position in general. The position de- 
pendence, however, is negligible for frequencies 
that are low in comparison with the characteristic 
frequencies of wave propagation in the chamber. 
As previously stated, these low frequencies are 
the onfy ones considered in the later sections of 
the report. 

In general, a precise analytical evaluation of 
the constants in equation (10) does not appear 
possible at the present time; the knowledge avail- 
able concerning the atomization and combustion 
processes is too limited. Some insight, however, 
regarding approximate values can be gained for 
particular cases. For vaporization-limited com- 
bustion, Spalding’s (ref. 17) modification of the 
Frossling evaporation equation (ref. IS) offers 
approximate values of m and n. The modified 
equation can be written as 

<1ZP 12 lii (1 + 7?) D(1+0 27c Sc m Re m ) 

(u pi oc 


from winch the exponents for the relation 


are found to he 


dr 


D m p n 


1+| 0.276 Sc 1 ' 3 Re l/ * 
m = 1+0.276 "S+ 3 /?e l/2 

1/ 0.276 Sc m W /2 \/ . AT T [U ±T/T\ 
H 2 Vl + 0.276 Sc l/3 rte l/2 ) V + A p/p A pjp ) 
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when the gas viscosity n, the transfer number B y 
and the Schmidt number Sc are treated as con- 
stants. Since the Schmidt number will be of the 
order of 0.7, it can readily be seen from the ex- 
pression for the exponent m that its value will be 
about 3/2 for high Reynolds number (Re ^$>16). 
Under like circumstances, the interaction index n 
will have a value near ){ provided the perturbations 
in relative drop velocity and gas temperature (or 
their net effect) are small in comparison to the 
perturbation in pressure. The evaporation equa- 
tion, of course, must be treated in more detail to 
obtain an approximate value of the steady-state 
time lag a. The more detailed approach was 
used in calculating the droplet histories presented 
in reference 19. These histories, which are given 
for several propellant combinations, offer a con- 
venient means of estimating the steady-state 
time lag. If the time lag is assumed to include 
most of the drop lifetime, the quantity y appearing 
in equation (10) will have a value near unity. 

The empirical equation of reference 13 for the 
mean drop diameter produced by two impinging 
jets can be used in evaluating the parameter a 
for injectors of this type. With the assumption 
of very low gas velocity in the region of impinge- 
ment, the equation to use is 

h/D — 2 . G4 \ hV 0.97 hV 

where* hole diameter h and drop diameter D are 
expressed in inches and injection velocity U is in 
feet per second. From this equation, the formula 
for a is found to be 

0.97 h,V+ 1 .32 -yJhV 

0.97 A F+2.64 \1,T 

Accordingly, the parameter a, which is termed the 
atomization index, will have values between % 
and unity, depending on the magnitude of the 
product hY. For the usual range of Jt\ (){ to 7), 
the formula indicates the atomization index will 
average about %. This is quickly ascertained 
from the plot of the formula in figure 1. 

If tin* atomization index a is put equal to zero, 
the term containing the perturbation in injection 
velocity disappears from equation (10), and the 
equation is reduced to the one used by previous 
authors. Such a procedure, of course, is tanta- 
mount to assuming initial drop size to be inde- 



h V, (in.)Cft /sec) 


FrorRK 1. Variation of atomization index a with product 
hV of hole diameter and injection velocity for two 
impinging jets (bast'd on empirical equation of ref. 13). 

pendent of injection velocity. In view of the 
observations of references 12, 13, and lb, however, 
this is not a generally valid assumption, and the 
complete equation must be used to define the rate 
of change of the time lag. It is t lie intent of t lie 
following analyses to show the important role the 
inject ion- velocity term has in the theory of cham- 
ber dynamics and combustion instability. In order 
to emphasize the relative effect of neglecting this 
term, (be analyses will be made with several 
simplifying st ipulat ions. 

CHAMBER EQUATION AND TRANSFER FUNCTIONS 

As brought out in the foregoing section, the 
time-lag concept provides a convenient combus- 
tion model consisting of a distribution of burnt-gas 
sources. As a consequence, the problem of 
chamber dynamics becomes one of a gas-flow 
system fed by these sources and depleted by the 
flow through the nozzle. In accordance with the 
principle of mass continuity, the excess of the 
feeding rate over the depletion rate must at each 
instant be equal to the rate of change of the mass 
within. In terms of small fractional perturba- 
tions, this can be expressed as follows: 


A uh 
w b 


, , A w e . , d 



(II) 


where the dimensionless time z is the ratio of the 
actual time to the average steady-state residence 
lime of the gas in the chamber. This gas 
residence time 6 g , of course, is equivalent to 
A I[w t and can be evaluated by multiplying the 
chamber volume by the factor pjwJVr, 

The perturbation &w b fw b in the feeding, or 
gas generation, rate represents the integrated 
effect of all the sources; the contribution of each 
source is dependent on the rate of change of its 
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associated time lag. The general problem in- 
volving nomniiform lime lags is rather involved 
alhegraically, but can be treated by using the 
methods of reference 11. For the present pur- 
pose, however, the time lag is assumed to be the 
same for all the sources. With this simplification, 
the perturbation in the total rate of gas generation 
can be related to the time lag in the following 
manner (ref. 11): 


^(2)=^ 

w b Wt d z 


( 12 ) 


where the time lag is expressed in the dimension- 
less form r, which is equivalent to ojB g . By 
using equation (10) and considering a bipropellant 
rocket in which the fuel drops have the longer 
lifetimes, the rate of change of the time lag is 
given by the equation 


dr /r> . — d AWf , —v 

— =— (3— m)yaT -r- — - {z~ t) 
d s ' dz w f 


— n 



(13) 


where, at each instant, the pressure is considered 
uniform throughout the chamber for the range of 
low frequencies being considered. The fractional 
change in the total injection rate can be separated 
into its oxidizer and fuel components by the 
relation 


Wi 


(2-t) 


r £w 0 
r + 1 w 0 


(* 


r)d 


1 Aw f 
r+1 w f 


(*-r) 

(14) 


Upon combining equations (12), (13), and (14), 
the equation for the perturbation in the rate of 
burnt-gas generation is found to be 


(*>-4- 4- ( 2 --r)+4r & (—’> 

w„ r 1 1 w„ r+1 Wf 


, _ d A w, , 

+ (3— m)yar ' (2—7) 


+„[4 w -f <*-?>] os) 


sumption can be made with reasonable confidence 
for many applications where the adiabatic flame 
temperature varies only slightly with mixture 
ratio in steady state. With this assumption, the 
dynamics of the unsteady nozzle flow can be 
accounted for in the low-frequency range simply 
by including a portion of the nozzle volume in the 
evaluation of the gas residence time (ref. 20). 
The nozzle flow, therefore, is treated as if it were 
steady. Finally, the specific-heat ratio is con- 
sidered to be unity. Accordingly, the perturba- 
tion in the flow of burnt gas from the combustion 
chamber is given by the equation 


A v> € 
w e 



(16) 


For a more detailed treatment of the exhaust 
problem, the reader is referred to the work of 
Cheng (ref. 20), Tsion (ref. 21), and Crocco (ref. 
22). (Crocco’s treatment is also presented as an 
appendix in ref. 11.) 

The simplification made in regard to the mass- 
perturbation term of the continuity statement is 
as follows. It is assumed that at each instant the 
integral over the chamber length of the entropy 
perturbation is nearly zero and can be neglected. 
The perturbation in the total mass of gas within 
the combustion chamber can be obtained, there- 
fore, from the equation 




(17) 


where the specific-heat ratio is again considered 
to be unity. A treatment retaining the integrated 
entropy change is not overly complicated if it is 
first assumed that the combustion is concentrated 
in a front near the injector face (ref. 11). Such a 
t reatmen t , however, exaggerat es the cut ropy 
effect if the actual combustion is spread at all in 
the axial direction. 

Upon combining equations (15), (16), and (17) 
with the continuity expression of equation (11), 
the combustion-chamber equation relating per- 
turbations in pressure and propellant flow is found 
to be 


In keeping with the intent of the present 
analysis, the flow leaving the combustion chamber 
is treated in the following simplified manner. 
First, it is assumed that the entropy perturbation 
at the nozzle entrance is negligible. This as- 


d^ !(z) + (1 _ B) fV (g)+B ^ (2 _ 7) ' ^ (2 _7) 

ds p w p P r- t-lw„ 

-7+1 U (s_T) + (3_m) yaT Tz W (z ~ 7 } 

(18) 
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Since the solutions of equation (18) are of the 
type e sz , the equation may also be written in 
operator form as 

r , , ,, w Ap r A w 0 - 

11+5— n(l — e“ rs ) -d-=— j-r -=+ e~ is 
1 p r+1 w 0 

++ [l + (r+l) +' <T 7 ’ (19) 

wliere the dimensionless operator s is, in general, 
complex. The factor rendering the operator 
dimensionless is, of course, the gas residence time. 

From equation (19), the following transfer 
functions of the combustion chamber can be 
written: 


1+5 — 71(1—6 T5 ) 


/ Ap/p_\ 

r 

\AW U /'(V U / c orn 

t .Uj r r i 6 

/ ±P/7> \ 

- 1 e -i 

\&WfJ st . 

«■„ 1 

ap!p \ 

_/ Ap/J) \ 

A1+/++ con at. r 

\At Vf/usf/t 


1+5— n(l — e Vs ) 


l + (3 —m)ycLTS , 

1+5— n(l— e~ fs ) ( } 


Despite the simplifying assumptions used in the 
derivation, it is believed these transfer functions 
describe the predominant form of response for 
many applications and thus are useful approxima- 
tions for controls work. It should be remembered, 
however, that in the present report the fuel drops 
are considered to have lifetimes longer than those 
of the oxidizer drops. If the situation is reversed, 
equations (20) and (21) must be revised as follows: 

r+1 

_ . H (3— m)yar* 

/ A pfp_\ _ r g _- rs r 

w 0 / cons i' Wf r+1 1+5 ft(l c *) 


(A 

\AwJwJ, 

/ Ap/£ 
\Mv f /w 


±P/P_ \ 

\W fj\V f / conbt . w Q 


1+5 — 71 ( 1—6 TS ) 


In the controls field, it is often convenient to 
consider dynamic relations in block diagram form; 
hence, a block diagram representation of the 
chamber dynamics is given in figure 2. Equation 
(19) is portrayed by the portion of the diagram 
drawn in solid lines. The term 

1 +5— n(l — 



Figure 2.-- Block diagram of chamber equation (propellant system feedback shown as dashed lines); 
534118—60 —2 
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of tin* equation is broken down into an equivalent 
closed loop in order to emphasize the positive 
feedback that exists within the combustion 
chamber. This feedback action, of course, results 
from tlu* sensitivity of the time lag to chamber 
pressure. Instability of this internal loop with 
constant injection rate is treated extensively in 
reference 11 under the designation “intrinsic 
instability.” In practice, of course, low-frequency 
combustion instability (chugging) also involves 
the feedback paths provided by the fuel and 
oxidizer systems. These external feedbacks, 
which are negative by virtue of the transfer 
functions involved, are indicated by the dashed 
lines in the diagram. 

FREQUENCY RESPONSE 

When the chamber t ransfer functions (eqs. (20), 
(21), and (22)) are normalized to unity gain, they 
(*an be considered in the general form 

1 + Qts 

£ T S 1 ^ 

1+5 — n (1— c~ r *) 

where the parameter Q is defined according to the 
following table: 


Normalized transfer function 

Q 

r r 1 / Ap/7) \ 

r \ A Wolfits const, wj 

0 

, , ,, / &PfP_\ 

\ AtCy/lt’// const. w Q 

(r+l) (3 - to) yet 

( AvFp \ . ( \ 

\ &W 0 jw 0 J const, r \ &WfjWf) const, r 

(3 — m) yet 


The parameter Q will be referred to as the 
“injection-velocity sensitivity” since its magnitude 
determines to what extent the response will be 
affected by the perturbation in injection velocity 
for a given value of time lag and interaction index. 

The frequency response defined by the general 
transfer function for various values of Q and r 
is presented in figures 3 and 4 at the end of this 
report. (The value of n is kept constant at 0.5 in 
these figures.) The amplitudes (fig. 3) and the 
phase shifts (fig. 4) are given as functions of the 
dimensionless frequency of neutral oscillation 0. 
This dimensionless frequency, of course, is equiva- 
lent to the dimensional angular frequency multi- 


plied by the gas residence time. (In using the 
computations presented herein, it is well to note 
that the imposed frequency limit can be expressed 
approximately as /3<1 /Mach number.) The range 
of values shown for the sensitivity parameter Q is 
believed to be a realistic one. Consider, for ex- 
ample, a combustion process conforming to the ap- 
proximate values previously discussed for the quan- 
tities m, y, and a (3/2, 1, and 2/3, respectively); 
the value of the sensitivity Q would then be unity 

ApIp._ 

M 

/ Ap/p \ 

\Aw f /w f J e 


for the transfer functions ( — 


s (±P%-) 
\Aw u/ !w 0 / c 


and 


and even larger for the transfer func- 


tion (r+l) (APJ”) 

/ conn. w 0 


All the transfer func- 


tions, of course, woul< 
for the zero-sensitivity case 


yield the 
r+ 1 


curves shown 

(Ap/p \ 

\ Aw? fin ) 
\^ w oi w o/ const, it j 

if the injection-velocity effect had been neglected 
in the time-lag concept. These curves, therefore, 
can be used as a basis of comparison for the more 
general time-lag concept, 

From this information presented in figures 3 and 
4, it can be concluded that the injection-velocity 
effect on the time lag is an important factor to 
consider when treating the dynamics of the com- 
bustion chamber. This is evident upon noting 
the strong dependence of the frequency response 
on the value of the injection -velocity sensitivity. 
At the higher values of the sensitivity parameter, 
both the amplitude and phase curves are con- 
siderably different than those for a sensitivity of 
zero. It should be noticed, too, that the difference 
becomes greater with increasing time lag. 


STABILITY CONSIDERATIONS 


For an analysis of combustion instability, the 
chamber equation (eq. (18)) must be considered 
together with appropriate equations for the fuel 
and oxidizer feeding systems. An equation for 
propellant -feed systems that can be applied to 
either pump- or tank-fed rockets is presented in 
reference 11. This equation, with certain simpli- 
fying assumptions, takes into account the resist- 
ance, inductance, and capacitance of the feeding 
system. In stability evaluations of particular 
systems, in which these effects are specified by 
given constants, the equation can be used with- 
out difficulty. A general treatment, however, 
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becomes very involved when the complete equa- 
tion is used. 

Therefore, in order to show without undue com- 
plication the relative importance of the injection- 
velocity effect, a simple, tank-fed system with 
negligible capacitance and inductance is assumed 
for both propellants in the presold treatment. 
Such a system represents the limiting case ap- 
proached by tank-fed rockets with short lines. 
The simplified equations for the fuel and oxidizer 
feeding systems can be written as follows (both 
tanks have the same pressure): 


Another equation is apparent after dividing the 
first of equations (26) by the second: 


,-1+^tan rfi R 

Vl + tftn 2 70 


(28) 


Introducing equation (28) into equation (27) 
(and rearranging) gives the following equation for 
7 in terms of p and the system parameters: 


- = ±\P 2j r ( n - 1 ) 2 — ( 71 + K) - 

QKp 


(20) 


the square root is negative if, and only 


where 


Av} i == _ k ±p 

(23) 

if, the parameter Q is negative. The second 

1C; p 

relation required between 7 and p is given by 

An Ap 

(24) 

equation (28), which can be put in the quadratic 
form 

n ~ v 




f (r? + K) 2 — p 2 } t an 2 7/3— 2 (n — 1 ) p t a n 7/3 

K= V 

2 U\.f-1>) 


+ (n+K)*-{n-\)*=0 (30) 


The el i min ant of equations (18), (23), and (24) 
yields the following characteristic equation: 


where 


l+8-n(l-e-*)+ Ke~' s ( 1 + QTs) - 0 (25) 

Q=(3—m)ya 


From equation (30), it is found that 
tan rp 


The characteristic equation determines the bound- 
ary between stable and unstable oscillation when 
the real part of the operator s is zero. The oper- 
ator, therefore, is set equal to ip, where p is the 
dimensionless frequency of neutral oscillation. 
By separating the real and imaginary parts of 
equation (25), two real equations are obtained: 


(n-1) P± (n + K) A p 2 + (n- 1 ) 2 — (n 1 E ) 2 
P 2 —(n + K) 2 

(31) 

where the square root again bears the sign 
of Q. Solving equation (31) for 7 gives the 
equation 


- 1 r 

r = P r 


] 


(32) 


QjKp sin rp n— 1 — (n+K) cos rP 
QrKp cos rp=—p-\-(n-\-K) sin 7 p 


(26) 


Equations (26) represent the stability boundary 
in parametric form. Because of the transcendental 
nature of these equations, a single equation of the 
stability boundary cannot be obtained. The para- 
metric equations, however, can be expressed in 
more convenient form. 

Upon squaring and adding equations (26), it 
can be ascertained that 

(QtKP) 2 =P 2 + (n- 1) 2 + (n+K)* 

-2(»+g)"~'+' al “ nW (27) 

Vi -Man 1 70 


r=7 1 (tTT— tan 1 

(n-\)P+(n + K) V p 2 +(n-\) 2 -(n + Ky 
P 2 - -(n + K)> 

where the principal value of the arctangent is used 

^<tan _1 <^ and a is an integer. For the 

present purpose, consideration is limited to the 
fundamental oscillatory mode, that is, the mode 
corresponding to the smallest time lag satisfying 
equation (32); thus, the integer a is given a value 
of unity. The stability boundary, therefore, is 
expressed by writing equations (29) and (32) in 
the following simultaneous form: 


-_± \P 2 + (ft — 1 ) 2 — (n 4- K) 


QKP 




r=d ir — tan 


(n— 1)< 3± (n 4- K) \& 2 + (n l) 2 - (« i A'! 2 ”! 

P 2 -(n-\-K) 2 J J 


( 33 ) 
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Equations (33) can be solved directly for two 
special eases. First, it is readily seen that, when 
the time lag is*constant (Q— 0^n~ 0) } the solution 
is as follows: 

T=p (ir-tarr 1 0) (34) 

where 

The second case is that of a time lag dependent on 
chamber pressure but independent of injection 
velocity (Q= 0, 0). The following [solution is 

then ol)tained: 

7= K’ r - ta „- , 1 f M ) (35) 

where 

V (h "h A ) ‘ — (ft — 1 ) * 

For the general case, equations (33) were solved 
numerically by iterative procedure with a high- 
speed digital computer. The stability boundaries 
obtained are presented in figure 5 at the end of 
this report. The critical values of the time lag 7, 
that is, the values giving neutral oscillation, are 
shown as functions of the injection-velocity sensi- 
tivity Q for various values of feed system pressure 
drop ( P 0t f—p)fp . This pressure-drop parameter, 
of course, is equivalent to the quantity 1/2 K. 
Separate charts are given for different values of 
the interaction index n. The unstable domain 
associated with each stability boundary is repre- 
sented by the region of the figure above the 
boundary. Curves of constant critical frequency 
£ (dashed lines) are also given in the charts. 

The conclusion drawn from figure 5 is that the 
injection-velocity sensitivity cannot be neglected 
in the theory of low-frequency combustion insta- 
bility. As shown in the figure, the critical time 
lag for a given interaction index and pressure drop 
is highly dependent on the magnitude of the 
injection-velocity sensitivity. The influence of the 
sensitivity parameter Q is shown to be generally 
destabilizing; however, it is also apparent that 
small values of the parameter Q can be stabilizing. 
Comparing the three parts of the figure shows that 
the destabilizing effect of the larger sensitivities 
is greater at the smaller values of interaction index. 

This latter observation is more evident from 
iigure 6, where, for a constant pressure drop, 



Figure C). Effect of interaction index n on stability 
boundary. Pressure drop (P 0 ,f — p)/p> 0.4. 


stability boundaries for different values of inter- 
action index are shown together. This figure also 
shows t lie influence of the interaction index to be 
small at the higher values of (he injection-velocity 
sensitivity. 

In the final figure (fig. 7), (lie critical curves 
for an interaction index of 0.2 are given for 



Figure 7. — Effect of injection- velocity sensitivity Q on 
unconditional stability. Interaction index », 0.2. 


constant values of injection- velocity sensitivity 
with the pressure drop and time lag as ordinate 
and abscissa, respectively. With this low value 
of interaction index, there is, according to reference 
11, a value of pressure drop large enough to make 
the S 3 7 stem unconditionally stable, that is, stable 
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for all values of time lag. The criterion of reference 
11 for unconditional stability can be written as 

(»<D 

For an interaction index of 0.2, the criterion yields 
a value of 0.83 for the pressure drop. This value 
is marked in figure 7 by the dashed line. Since 
the line is crossed by stability boundaries, it is 
evident that the criterion is not sufficient with the 
generalized time-lag concept. In fact, equation 
( 29 ) indicates that unconditional stability does not 
exist [for injection-velocity sensitivities other 
than zero. This can be shown upon writing the 
equation as follows: 


r, 2 f A T 1 ) (2/i T K 1 ) 

1 — (Qr AO 1 2 3 4 5 


( 36 ) 


Equation ( 36 ) determines the following two 
conditions having no real frequency of neutral 
oscillat ion : 


(1) K< 1 —2 n and K<C 7 r= 

Qt 



( 2 ) K> 1 


-2 n and K^> ’77= I 

M r \ 


The first condition is the criterion for stability 
(the second for instability). Tn terms of pressure 
drop, it is expressed as follows: 


r r<,'f-p 

p 

w v 




The criterion, although sufficient for stability, 
does involve the time lag; the stability, therefore, 
is not unconditional. 

CONCLUDING REMARKS 

In conclusion, it is important to recall that the 
time-lag concept used heroin was introduced 
because previous studies indicated a dependency 
of combustion time lag on injection velocity. 
The purpose of the present study was to show the 
importance of using this more general concept 
when treating chamber dynamics and stability. 

The analyses made with the generalized concept 
of time lag have shown that chamber responses 
and system stability boundaries can be greatly 
affected in the low-frequency range by reasonable 
variations of the time lag with injection velocity. 
It can be concluded, therefore, that the injection- 
velocity effect on the time lag is an important 
factor to consider in the theory of chamber 
dynamics and combustion instability. 

In view of this conclusion, it can be stated that 
information concerning the relation of drop size 
to injection velocity is essential for the accurate 
prediction of rocket dynamics and stability. 
Investigations, therefore, should be conducted 
to provide this information for the various injector 
types of interest. 


Lkwir Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, April 24 , 1059 
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Dimension! ess frequency, ft 


(b) Dimensionless time lag r, 0.6. 

Figure 3. Continued. Amplitude of frequency response for general form of normalized transfer functions. Interaction 

index n t 0.5. 









Amplitude 



(c) Dimensionless time tag t, 0.8. 

Figure 3. — Continued. Amplitude of frequency response for general form of normalized transfer functions. Interaction 

index n, 0.5. 
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(d) Dimensionless time lag r, 1.0. 

Figure 3. Continued. Amplitude of frequency response for general form of normalized transfer functions. Interaction 

index u f 0.5. 
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Dimensionless frequency, /? 


(c) Dimensionless time lug r , 1.2. 

Figure 3.- Continued. Amplitude of frequency response for general form of normalized transfer functions. Interaction 

index n, 0.5. 
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Figure 3.- Concluded. Amplitude of frequency response for general form of normalized transfer functions. Interaction 

index n, 0.5. 
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Dimensionless frequency, /3 

(d) Dimensionless time lag 7, 1.0. 

Figure 4. — Continued. Phase shift of frequency response for general form of normalized transfer functions. Interaction index n } 0.5 
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(f) Dimensionless time lag r, 1.4. 

Figure 4, — Concluded. Phase shift of frequency response for general form of normalized transfer functions. Interaction index n, 0.5. 
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(a) Interaction index to, 0.5. 

Figure 5. : — Stability boundaries of tank-fed rocket with negligible feeding system inductance and capacitance. 







Critical value at time lag 
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(b) Interaction index n, 0.3. 

Figure 5. — Continued. Stability boundaries of tank-fed rocket with negligible feeding system inductance and capacitance, 
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Figure 5. — Concluded. 


(c) Interaction index n, OJ. 

Stability boundaries of tank-fed rocket with negligible feeding system inductance and capacitance. 
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